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In this work we consider the universal crossover behavior of two non-equilibrium systems exhibiting 
a continuous phase transition. Focusing on the field driven crossover from mean-field to non-mean- 
field scaling behavior we show that the well-known Widom scaling law is violated for the effective 
exponents in the so-called crossover regime. 
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I. INTRODUCTION 

A foundation for the understanding of critical phe- 
nomena was provided by Wilson's renormalization group 
(RG) approach 0, which maps the critical point onto 
a fixed point of a certain transformation of the sys- 
tem's Hamiltonian, Langevin equation, etc. The RG 
theory presents a powerful tool to estimate the val- 
ues of the critical exponents, it allows to predict which 
parameters determine the universality class and it ex- 
plains the existence of an upper critical dimension Dc 
above which the mean-field theory applies. Furthermore, 
crossover phenomena between two different universality 
classes are well understood in terms of competing fix- 
points (see for instance f3l)- Nevertheless there are still 
some open aspects of crossover phenomena which are 
discussed in the literature. The question whether the 
crossover scaling functions are universal was revisited 
several times 1^ |5|, ,7,, ,^ 0, 0| . For instance it was 
shown recentlty that two different models, belonging to 
the same universality class, are characterized by the same 
(universal) crossover scaling functions [Tol |. This is re- 
markable, since the universal scaling behavior is usually 
restricted to a small vicinity around the critical point. 
In the case of a crossover the universal scaling functions 
span several decades in temperature or conjugated field. 

Another question of interest concerns the so-called ef- 
fective exponents jlll | which can be defined as logarith- 
mic derivatives of the corresponding scaling functions. 
It is still open whether these effective exponents fulfill 
over the full crossover the well-known scaling laws which 
connect critical exponents. This question is closely re- 
lated to the more general and very important question 
whether effective exponents obey the scaling laws at all. 
For instance it is known experimentally |l2j | as well as 
theoretically 0| that the asymptotic scaling behavior 
is often masked by corrections to scaling, so-called con- 
fluent singularities. In this case it is useful to analyze 
the data in terms of effective exponents and the above 
question naturally arises Thus the validity of the 

scaling laws for effective exponents was addressed in ex- 
perimental works, RG approaches as well as numerical 
simulations. In particular the violation of the scaling 
laws for effective exponents was conjectured from a RG 



approximation 15]. But neither experimental nor nu- 
merical work could clearly confirm this conjecture so far. 
For instance Binder and Luijten considered numerically a 
crossover in the Ising model and discussed the validity of 
the Rushbrook scaling law |^. The observed nonmono- 
tonic crossover behavior suggests again a violation of the 
Rushbrook scaling law. However it can not be considered 
as a rigorous proof. 

In this work we consider a field driven crossover (in 
the so-called critical crossover limit 0, 0) in two dif- 
ferent models exhibiting a non-equilibrium second order 
phase transition. Varying the range of interactions we 
investigate the crossover from mean-field to non-mean- 
field scaling behavior. The order parameter and the or- 
der parameter susceptibility is measured as a function of 
the conjugated field and we are able to determine the 
corresponding effective exponents over the full crossover 
regime. Our results show that the well-known Widom 
scaling law is clearly violated for the effective crossover 
exponents. Furthermore we present a simple analytical 
argument, suggesting that the scaling laws are valid for 
the asymptotic scaling regimes (where the systems are 
characterized by a pure algebraic behavior) , whereas the 
scaling laws do not hold for the crossover regime (char- 
acterized by a non algebraic behavior). 

II. MODELS AND SIMULATIONS 

In the following we consider two different cellular au- 
tomata exhibiting a so-called absorbing phase transition. 
The first model is the conserved transfer threshold pro- 
cess (CTTP) 0. In this model lattice sites may be 
empty, occupied by one particle, or occupied by two par- 
ticles. Empty and single occupied sites are considered 
as inactive whereas double occupied lattice sites are con- 
sidered as active. In the latter case one tries to transfer 
both particles of a given active site to randomly chosen 
empty or single occupied nearest neighbor sites. 

The second model is a modified version of the Manna 
sandpile model the fixed-energy Manna model [l8l| . 
In contrast to the CTTP the Manna model allows un- 
limited particle occupation of lattice sites. Lattice sites 
which are occupied by at least two particles are consid- 
ered as active and all particles are moved to the neigh- 
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boring sites selected at random. 

In our simulations (see p^ . [20j for details) we have used 
square lattices (with periodic boundaries) of linear size 
L < 2048. All simulations start from a random distribu- 
tion of particles. After a transient regime both models 
reach a steady state characterized by the density of active 
sites The density is the order parameter and the 
particle density p is the control parameter of the absorb- 
ing phase transition, i.e., the order parameter vanishes at 
the critical density pc according to, p^ oc 5p^ , with the 
reduced control parameter dp — p/ pc^ 1. Below the crit- 
ical density (in the absorbing phase) the order parameter 
is zero in the steady state. 

Similar to equilibrium phase transitions it is possible 
in the case of absorbing phase transitions to apply an ex- 
ternal field h which is conjugated to the order parameter. 
The conjugated field has to act as a spontaneous creation 
of active particles, destroying the absorbing state and 
therefore the phase transition itself. Furthermore, the 
associated linear response function Xa = dp^/dh has to 
diverge at the critical point {5p = 0, h = 0). A realization 
of the external field for absorbing phase transitions with 
a conserved field was recently developed in where 
the external field triggers movements of inactive parti- 
cles which may be activated in this way. At the critical 
density p^ the order parameter scales as p^ oc h^^°' . Us- 
ing the conjugated field it is possible to investigate the 
equation of state p^{p,h), i.e., the order parameter as a 
function of both the control parameter and the external 
field. A recently performed scaling analysis reveals that 
the CTTP and the Manna model are characterized by the 
same critical exponents as well as by the same universal 
scaling form of the equation of state, i.e., both models 
belong to the same universality class |2ll |. 

According to the above definition particles of active 
sites are moved to nearest neighbors only, i.e., the range 
of interactions is i? = 1 . It is straightforward to imple- 
ment various ranges of interactions into these models [10| . 
In this modified models particles of active sites are moved 
(according to the rules of each model) to randomly se- 
lected sites within a radius R. Since the dynamics of 
both considered models is characterized by simple parti- 
cle hopping processes, various interaction ranges can be 
easily implemented and high accurate data are available. 
This is a significant advantage compared to e.g. equi- 
librium system like the Ising model where the increasing 
interaction range causes a slowing down of the dynamics. 

For any finite interaction range the phase transition is 
characterized by non-mean-field scaling behavior which 
now takes place at the critical density p^ rt- A mean-field 
phase transition occur for infinite interactions {R oo) 
only. But mean-field behavior could occur away from 
the critical point if the long range interactions reduce 
the critical fluctuations sufficiently. The crossover be- 
tween the mean-field and non-mean-field scaling regimes 
is described by the famous Ginzburg criterion 22] which 
states that the mean-field picture is self-consistent in the 
active phase as long as the fiuctuations within a correla- 



tion volume are small compared to the order parameter 
itself. This leads for zero field to the crossover condi- 
tion 0[R^ii [p — p r )*^] = 1, with the crossover exponent 
= (4 — D)/2D ITo]. In order to avoid lattice effects we 
use the effective range of interactions 23j 

Rl = -Efc. \L.-rJ<R. (1) 

where z denotes the number of lattice sites within a ra- 
dius R. 



III. UNIVERSAL CROSSOVER SCALING 

The crossover scaling function has to incorporate three 
scaling fields (the control parameter, the external field, 
and the range of interactions), i.e., we make the phe- 
nomenological ansatz 



p,{p,h,R,ft) ^ 



(2) 



where the universal scaling function $H is the same for all 
models belonging to a given universality class whereas all 
non-universal system-dependent features (e.g. the lattice 
structure, the update scheme, etc.) are contained in the 
so-called non-universal metric factors Op, a^, and [23| . 
These factors are determined by three conditions which 
normalize the scaling function 91. Firstjthe analytically 
known mean-field scaling function |2lll25l| 



RMF{x,y) = - 



(3) 



should be recovered for R ^ oo, i.e., 9i{x,y,0) = 
RMAx,y)- Therefore, 9\{1,0,0) = ^mf(1,0) = 1, 
5H(0, 1,0) = i?MF(0, 1) = 1, implying = Up^R^^/ p^^n^^ 
and Oh — a,i,j,^3o. Finally, the non-universal metric fac- 
tor Or can be determined by the condition [H(a;, 0, 1) ~ 
for X ^ Q yielding [lo| 



Pc.R=l ^p,i?.^c 



3d/(/3mf-/3d) 



(4) 



The metric factors were already determined in previous 
works 0, I23, thus no parameter fitting is needed in 
order to perform the following scaling analysis. 

In this work we focus our attention to the field driven 
crossover, i.e., we consider the CTTP and the Manna 
model at the critical densities p^ ^ which were determined 
in 10]. In Fig. n we plot the corresponding data of the 
CTTP for various values of the interaction range R. As 
one can see the power law behavior of the order param- 
eter changes with increasing range of interactions. 

The scaling form at the critical point is given by (set- 
ting a-^RZ-} = 1) 



p^{p^,R,h,R,tf) - 

(aRi?.„.)"''"^^'^ 9^(0, a,/i aT^''^R:r'^, 



(5) 
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FIG. 1: The universal crossover scaling function 9^(0, a;, 1) 
of the order parameter of the CTTP and Manna model at 
the critical density for D = 2. The metric factors are 
given by e = ana'^ and d = a^. The dashed lines cor- 
respond to the asymptotic behavior of the two-dimensional 
system {l3o=2/ 0-0=2 = 0.287) and of the mean-field behav- 
ior (/3mf/o"mf = 1/2). The universal amplitude is given by 
771a, h = 0.681. The inset displays the order parameter of 
the CTTP for various values of the interaction range R £ 
{1, 2, 4, .., 128} (from top to bottom). The dashed lines are 
just to guide the eyes. 



with /3mf = 1 and CTmf — 2, respectively. For sufRciently 
small field the universal function scales as 



5H(0,a;,l) 



for 



0. 



(6) 



with the universal amplitude m^^h- The scaling form 
Eq. jSJ has to equal for R=l the D-dimensional scaling 

ansatz [ ^ (a^^^i/i)''"/'^" ] leading to 



/5mf/0— o"mf/3d /cd 



(7) 



According to the scaling form Eq. © we plot in Fig.Q] 
the rescaled order parameter pa (oni?„ff)^ as a function of 
the rescaled field ah/i(aRi?„ff)^. We observe an excellent 
data collapse for the full crossover behavior confirming 
the above phenomenological scaling ansatz. 

However, since the entire crossover region covered sev- 
eral decades it could be difficult to observe small but 
systematic differences between the scaling functions of 
both models. Therefore, it its customary to scrutinize the 
crossover via the so-called effective exponent [gII^ IToIITH 



d 



dlnx 



ln5n(0,x, 1). 



(8) 



The perfect collapse of the corresponding data is shown in 
Fig. Inland confirms again the universality of the crossover 
scaling function $H. 
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FIG. 2: The effective exponent (/3/a)eft. 

Next we consider the order parameter susceptibility. 
The scaling form of the susceptibility is given by 

a^Xa(/0, ^ (9) 

AT"- €{a,{p - p.,,) A, a,h A^-^ a-ii?"/ A^) . 

On the other hand the susceptibility is defined as the 
derivative of the order parameter with respect to the con- 
jugated field 



d 

X.{p,h,R^!!) = —p^{Sp,h) 



(10) 



~ a,A^"--'^«-5n'(a,(p - p.,«)A, a,/^A^«^ a-ii?-iA*) 

with 9l'(a;, y, z) — dyyi{x,y, z). By comparing this ex- 
pression with Eq. we find €{x,y,z) = dyd\{x,y, z), 
— a^^, as well as the Widom scaling law 



a - (3 



ill) 



which is well known from equilibrium phase transitions. 

Again, the mean-field behavior is recovered for R 
«3, i.e., £(x,2/,0) = C»^{x,y) = 1/2 {y + {x/2f)-^/^ , 
implying £(1,0,0) = Cmf(1,0) = 1, £(0,1,0) = 
CmfIO, 1) = 1/2, as well as 7mf — 1. 

Similar to the order parameter we plot the susceptibil- 
ity according to the scaling form 



(aHi?„„)T"-/^ £(0, a^h (a^R^,,) 



(12) 



Approaching the transition point the susceptibility is ex- 
pected to scales as £(0, x, 1) ^ m^ ,,x^'^"^'^" , for x ^ , 
where the universal power-law amplitude is given by 



'■h,R=i \ — 7MF/</'+crMF7r>/o"r>0 f-^o 

' I Wis 



(13) 
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FIG. 3: The universal crossover scaling function £(0, x, 1) 
of the susceptibility of the CTTP and the Manna model 
at the critical density for D — 2. The metric factors are 
given by e = and^ and b — a^^a|. The dashed lines cor- 
respond to the asymptotic behavior of the two-dimensional 
system (^'^d=iI o o=i = 0.713) and of the mean-field behav- 
ior (7mf/o"mp = 1/2). The universal amplitude is given by 
m^,h = 0.208. The inset displays the corresponding effective 
exponent {'^ja^^n. 



FIG. 4: The violation of Widom scaling law [Eq. 115|l 1 in the 
crossover regime. 



Using 1 — dXvLaxjdXvix we get 



ln9,5n(0,a;,l) = In ax - ln$K(0,a;,l) + c, (17) 



The rescaled susceptibility is shovi^n in Fig.O Over the 
entire crossover region we got an excellent data collapse 
including both asymptotic scaling regimes. The inset dis- 
plays the effective exponent 

(-) = -|^ln£(0,x,l) (14) 
V (T / ctf a In x 

which exhibits again a monotonic crossover from the two- 
dimensional scaling regime to the mean-field scaling be- 
havior. 



where c is some constant. It is straightforward to show 
that this differential equation is solved by simple power- 
laws [9^(0, a:,l) — CqX'^i with Ci — 1/aexpc]. Thus the 
Widom scaling law is fulfilled in the asymptotic regimes 
only. In the case that the scaling behavior is affected 
by crossovers, confiuent singularities, etc. no pure power- 
laws occur and the scaling laws do not hold for the cor- 
responding effective exponents. 



IV. WIDOM SCALING LAW 

In this way we have obtained the effective exponents 
(/3/(T)eff and (7/cr)ctf for the field driven crossover from 
mean- field to non- mean-field behavior. Thus we are able 
to check the corresponding Widom scaling law 

(-) -l-(-) ' (1^) 

for the whole crossover region. The corresponding data 
are shown in Fig. 21 As can be seen the Widom scal- 
ing law is fulfilled for the asymptotic regimes {D = 2 
scaling behavior and mean-field scaling) but it is clearly 
violated for the intermediate crossover region. This re- 
sult is not surprising if one notices that the above Widom 
law [Eq. H15|) ] corresponds to the differential equation [see 
Eqs. (IHlini)] 

- Ir^ IT ^(0, a;, 1) - 1 - ^^n(0, X, 1) . (16) 
avo-x ox omx 



V. CONCLUSION 

In conclusion, the crossover from mean-field to non- 
mean-field scaling behavior is numerically investigated 
for two different models exhibiting a second order phase 
transition. Increasing the range of interactions we are 
able to cover the full crossover region which spans sev- 
eral decades of the conjugated field. The corresponding 
data show that the Widom scaling law is violated in the 
crossover regime. Notice that we focus in our investi- 
gations on the particular universality class of absorbing 
phase transitions only for technical reasons. The demon- 
strated violation of the Widom scaling can be applied to 
continuous phase transitions in general. 
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